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DESIGN OP A COMPUTER AIDED OPTIMAL SPEED REGULATING SYSTEM 
FOE A BAR MILL MAIN DRIVE 

In this thesis an optimal design of a speed regulating system for 
a bar mill main drive has been obtained. 

In earlier designs of speed regulating systems, although concepts 
of mathematical modelling -were used, the merit of the control 
configuration ms judged by experimental observations and analog 
computer simulation. This left open the question of optimality. 

In this thesis, a new technique of designing speed regulating 
systems is evolved by using dynamic programming concepts. Use 
of dynamic programming enables one to answer questions about the 
goodness of design with certainty. 

There are two notable features of modelling of a main drive. One 
is that, for both armature current and air gap flux varying with 
time, a product term appears in the state equations. Another is 
that, the field system in main drive has a nonlinear B-H 
characteristic. 



It wuld appear that no explicit optimisation technique has 
so far been suggested for systems described by state equations 
that involve product terms. However, if one of the two variables 
(armature current and air gap flux) is time invariant, then a 
single valued piecewise linear approximation of B-H curve may be 
used. Consequently, optimisation techniques for "nearly linear" 
systems may be employed in this case also. 

The optimisation method developed here is based upon the premise 
that the product term can be approximated by a weighted linear 
combination of both the terms with time varying coefficients. 
Starting with an initial guess on the trajectories of both the 
variables and using them as time varying coefficients, the corres- 
ponding set of time varying linear or non-linear differential 
equations is constructed. The approximation of product terms is 
improved iteratively by taking Gateaux gradient functional of a 
non-negative function along each of the guessed trajectories. 

The optimal solution is obtained by using parametric expansion 
in each iteration. Thus one essentially solves a sequence of 
optimal problems vhich converge to the required solution. 

A control program implementing this optimisation process has been 
written in a form that is adaptable for use in control computers. 
On-line execution of this program for both generator field control 
mode as well as motor field control mode of operation has also 


been discussed 



CHAPTER - I 


INTRODUCTION 

1 2 

Following the publication of works of Pontryagin and Bellman , the 
term "optimum-design" acquired definite meaning in theoretical design 
of control systems. The engineering translations of these designs 
seem nearly feasible*, only when a fast interacting computer system 
is used as a control tool. 

This thesis concerns with design of a computer-aided-optimal speed 
regulating system of a certain bar mill drive. The optimality of 
operation is ensured by using Bellman's dynamic programming concepts in 
evolving the control logic. It is assumed that for on-line operations 
a process control computer would be available. 

In what follows an over-view of the thesis is provided. 

12 3 

1.1 Logic for System Design * 

The logic often used in attempting optimum system design iss 
if feasible solutions exist (satisfying design specifications) then 
there exists an optimum solution which not only satisfies the design 
specifications but is also optimiun with respect to some criterion. 

This argument provides the basis for the present investigation. 


* nearly ... because of finite measurement delays and processing 
times involved. 
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1 ,2 Feasible Solutions 


Both a.c, and d.c. drives have been employed in steel rolling mills. 
However, experience ’ has proven that the d.c. drive systems are to 
be preferred. Among the d.c. drive systems there are two arrangements: 
one has a three machine complex -vdiile the other one has one nachine and 
heavy rectifying eq^dpment accompanying it. The present investigation 
concerns a three raachine configuration as depicted in Pig.l-''. 


The two modes of operation^’ for the arrangement in Fig.lrl are 

1 ) Below the base speed of operation the error signal (difference 

between desired and actual speeds) is fed to the channel 
corresponding to generator field control, 

2 ) Above the base speed the error signal is fed to the channel 

corresponding to motor field control. 


The magnitude of the terminal voltage is utilised to exercise 
switching of error signal between the two modes of control. All the 
schemes, cited above, confirm to this arrangement. 

In a recent study”^ (Appendix A) on a control scheme similar to 

6 

Haskett's , the following was observed: 

a ) On employing field forcing units which are capable of supplying 
high d.c. voltage for short periods, one can obtain fast field 
current response. This can be explained by observing the 
solution of a first order differential equation as detailed in 
Pig.1-2. 




basic control configuration of three machine complex. 




00 0,5 1.0 1.5 2.0 

time in seconds ^ 

1) V(t) = u(t) 

2) Y(t) = 10 ^ U(t) t < 0.5 , V(t) = Tl(t) for t >0,5 

The model of field coil t- . 



Figure 1-2 

Current response of a field coil mth and without Tsormal excitation. 
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b ) It is easy to achieve the separation of two 


j 3 od.es of control 


'^7 appropriate hard?j-are. 


The studies cited above demonstrate that feasible 
the question v/hether these are also optimal^ in 


solntions exist bni't:;; 
, oiue s ens e , remains 


un---nswered* 


1 . 3 Investigation for Optir-Tuiii 

y V -i bar mill main (iir ^^•-itor0 

Having confirmed (App>endix A.; that for a part i on 

- ^r^stigate if it was 

feasible solutions existed, it was decided to 

. . control system.. 


I)ossible to design a computer aided optimal spe^^ 
uso of computer as a control aid is suggested 
\sTuld be available, as such, for economic reason^ 


stXLse a computer system 
9 , "1 0 


I?he Modelling of Mill Main Brive 


The mill drive has an overall dynamics of four'fcl’- 
in different time slices, either generator fiel^ 
field current is constant, one essentially has ^ 
The state equations of the process, when genera-'^*- 


ox’der, but because 


ovirrent or motor 


-t^ird order plant . 


(^xr control is 


-fcli® generator field 


exercised, are linear if one agrees to linearis^ 

, V nrtxe state equations 

B-H characteristics (a reasonable assumption^ 


of the process, when motor field control is exei'^^‘ 
and involve both a product type of non-1 inear i'ty 
functional non-linearity of the B-H curve. As 


j.sed, are non-line 


g,s well as the 

, field weakening/" 


X 4.V, • / XU V 'i • • S' 

strengthening (as the need be J is exercised, 

, ^ X ^ ^c:ti-onal non-linearity 

(usual range being O.5 p.u. to I.5 p.u.j, the 

of motor field systems cannot be over siinplifie*^-- 



PerfornancG Index 


For either mode of operation, the criterion of performance was chosen 
of the form 


P.I, 


0> L- 


1 A * (N - F + B * +-v/ + D * (V^ - 1.5.1 


.'t J-l T 

m md 


td' 


Fle-xihility in selection of the coefficients A, B and D is to he 
ohserved'*^''*^^"*^’'*^ to finally yield an acceptable design.. 


Opt imisation Procedure 

For the generator field control (linear control problem), Merriam's"**^ 
pa-rametric expansion technique has been used. 

To compute optimal control law for the motor field control mode of 
operation a new technique has been suggested in this thesis . The 
method requires replacement of product terms in the state equations 
by a. weij^ted linear combination of both the terms as explained in 
equation 1.5.2. 

x^(t) ^ X2 (t) = w * x^g('t) * X2(t) + (1 - w) * x^(t) * 
where 0 and ) 1.5.2 

x^g, X2g refer to guesses to be improved iteratively ) 

Now the optimal solution of the system equations with trajectory 
approximation as in equation 1,5.2, are obtained using parametric 
expansion method. Based on the solution the approximation in 1.5.2 is 
modified. The trajectories corresponding to (^ig» ^2g' “*’* are 

nanipulated by using Gateaux derivative of a positive-semidefinite' 




functional for each t in the interval rOjT]. The Gateaux differential 
of this fimctional equals zero over the entire time interval tdien the 
approximation is exact. One major advantage of this approximation is 
that Merriam's parameteric expansion can be used to obtain two initial 
value problems from the corresponding Hamilton Jacobi equations and 
solution of a two point boundary value problem is avoided. 

1,4 Main Contributions of the Thesis 

The object set for the thesis was to evolve an optimal design for a, 
bar mill main drive speed control system and this has been achieved. 

The optimal solution is obtained for continuous generator/motor field 
voltages, as it was assumed that continuous voltages would be available 
as forcing functions. 

In the course of this investigation a powerful trajectory manipulation 
technique, which can be easily implemented, using a digital computer, 
has been developed. The technique may be employed to optimise control 
problems characterised by multiplicative and functional (single-valued) 
non-linearities in system equations. Adequate proofs for the algorithm 
on optimising procedure are also provided. 

'*•5 Outline of other Chapters 

In Chapter II, the mathematical model associated with the bar mill 
main drive has been derived. The performance requirements on speed 
control system are spelled out next and a complete problem statement 


follows immediately 



To frame an appropriate context for the method suggested for the 
problems stated in Chapter II, a brief review of many algorithms 
available in the literature is presented in Chapter III. In the 
last two sections of the chapter, the new trajectory manipulation 
technique is discussed with complete details. Definitions and 
theorems concerning the proposed technique are also given. 

Chapter IV primarily concerns with the computational aspects of -fhe 
optimisation technique suggested in Chapter III for the problem 
stated in Chapter II. All the flow diagrams corresponding to 
computer program are described in detail. The computer results 
detailing the performance of the mill main drive are also given. 

The chapter ends with a discussion on on-line execution in a 
hypothetical situation, where a process control computer is 
available as a control tool. 

In the last chapter the main conclusions from this investigation 
are highlighted and directions for further research are suggested. 

The Appendix A, B and C contain material which is perljheral to 
main stream of thought in the thesis.. 



CHAPTER - II 


THE PROBLEM STATEMENT 

The purpose of this chapter is to obtain a mathematical model of the 
bar mill main drive and make a statement about the desired performance. 

f 

The modelling procedure employs an oriented system graph for the imin 
drive. The cutset and circuit equations for the graph are substituted 
in the terminal equations to finally yield the state equations. 

The design requirements are based on an industrial project report^*^. 
These specifications essentially define an admissible design. 

To achieve an optimum admissible design, a criterion of performance 
is chosen and a logic of iteration is followed to arrive at appropriate 
weighting factors to bo associated with each of the constituents of the 
performance criterion. In the same section the protection scheme for 
the main drive is also discussed because the logic of iteration 
(at computational stage) is ba&©dr-en protective features (at operational 
stage) tr.' be finally associated with the mill main drive. 

2,1 Derivation of State Model .. 

The state model for the system can be obtained by following graph- 

52 

theoretic procediire suggested by P.H.O'N.Roe . The method involves 
systematic substitution of circuit and cutset equations in the 
terminal equations to yield state equations for the given system. Bie 
forest is chosen so that most voltage drivers and capacitors are 



included in the forest while nost of the ciirrent drivers and 
inductors are consigned to cc-forest. For non-degenerate circuits, 
capacitor voltages and inductor currents finally turn out to he 
state variahles and, therefore, should he preferably measurable. 

In what follows the chosen forest is referred to as the formulation 
forest, the tree branch voltages and chord currents are called 
primary variables and the rest are called secondary variables. 

The system graph corresponding to a mill drive is shown in Fig. 2-1. 
The formulation forest is detailed by the solid lines in Pig. 2-1. 
The circuit and cutset equations can he written by arranging the 
primary and secondary variables to yield equations 2.1*2 having 
the form 2.1.1. 

Circuit Equations 


Cutset Sq-uations 


where A and T refer to across and through variables and 

subscripts p , s are used to distinguish primary/secondary 




A 


= , 0 


T 


= 0 


2 . 1 . 


variables. 
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Novj- on writing the terminal equations explicitely in primary 



Using circxiit and cutset equations we can eliminate secondary 
variables from the terminal equations 2.1.5 so that 



2 . 1 .5 









ir-7 


.!^2, 


Ki2 *^(13) 

. _ _■ 


K, 5 *'^(14) *66 

^14 

! 

1 

-Kij 


rl I^/dt 




y - Y 

18 

d I^/dt 
d I^/dt 




: / ^ 5 ! 


7 - V 

2 ^ _ 11 ^ 

vj-'y-'vyv^o’-'y' 

d Og/dt 



] 

-^4 * ’'15 ^ 


2 . 1.4 


From the equations 2 . 1 - 4 , the state equations (equation 2 . 1 . 5) can 
he derived 

d Ij/dt = - Tg * j 

d I^/dt = ^ / ^4 ^4 I 

d I^/dt = -(r^ + t^q) * / (L^ + L^) “ * ®6 * ^^^4^ ^ ^7^ ) ^'"*•5 

+ K^2 * ^ 

d §g/dt = K^5 * /(I4) *15/^6 - ®15 * ®6 / “"e - ^16 / 6 I 

2.2 The Uomalised State Equations 

The following system data is assumed to he given s 

G-en. ratings i 2040 -kw, 400-volts, J^O rpm 

r = r = O.OOI65 ohms, = Lc = O.IO75 mh 
ag 9 ®€ 2 

“ ^8 “ "'*26 ohms, = 4.4 h 

rated field current 50 amps. 










Motor ratings s 1875 kw, 4OO volts, 8O/16O rpm 


= ^10 0-00299 ohms, = O.315 mh 

“ "'ll ■ ''•° ^fm ' h ■ 

rated field current 100 amps 

2 

rotor inertia 10,000 kg. m 

2 

mill inertia 2820 kg.m 

2 

Total mill inertia = Jg = 12820 kg.m 
Prom the given system data normalised state equations can he obtained 
by defining following per unit* quantities. 

For the generator field circuit s 

Per unit voltage = 47»5 volts, Per unit current = 5O.O amps. 
Time constant = (l/0.288) seconds 

For the armature circuit s 

Per unit voltage = 386 volts, Per unit current » 5'I40 amps. 
Time constant = (1/1I.O9) seconds 

For motor field circuit s 

Per unit voltage = 100 V, Per unit current = 100 amps. 

Time constant = (1/0.2I8) seconds 

For the mechanical circuit J 

Per unit speed = 80 rpm 

Per unit torque = p.u. power/(p.u- rps x 2 ) 

= 256096 newton-meters 

Per znait inertia = p.u. torque/p. u. acceleration 
* 28O95.42 kg.m^ 

Total mill inertia is then 0,446 p.u. 

* ELseiidiere in the thesis per unit quantities are also referred to 
as base quantities. ■ 
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And the normalised state equations can be written as ■ 


d I^g/dt = - 0.288 * I^_ + 0,288 * V, 


■fg 


d I^/dt - 0.218 ♦ 


fg 

H. 0.218 . 


d I./dt . - 11.09 » I + 11.09 ♦ I, - 11.09 ♦ ^(I, ) • H„ 

^ 3. Iff m 

X 

d N/dt = 2.24 * - 2.24 * '^j, 

K- = I^/ ^(I^„) 


) 


2 . 2.1 


An asstunption of single valued relationship between (jZL , I. ) and 

ig fg 

> d ) leads to the choice of I- and I. as measurable state 
fm fm fg fm 

variables besides simplifying the overall structure of state 

equations*. In equation 2.2.1 the characteristic is 

assumed .to be linear while the , I„ ) characteristic is assumed 

fm fm 

to be piece-wise linear a.djusts fid. time const, in each segment). 


For either generator field control mode of operation or motor field 
control mode of operation the state equations can be obtained from 
equations 2.2.1 by observing that 

1 ) For generator field control mode 

^ ^ ^ ** '■ ** " 2 . 2.2 

2 ) For motor field control mode 


d I^g/dt » 0 .. .. .. .. 2.2.3 


* The double valued approximation of l) characteristics require that 
fluix be considered as a function of both I and (d l/dt). A hyateas^sis 
curve can be simulated by integrating an equation of the following form s 

d jZ!/dI ^ 1 / (1 - a * + b * I^ ) + c * d 3/d-fc, 

where a, b and c are constants. Constants a, b and c can be obtained 
from J2f(l, t) at t = 0, l(o) and d l(0+)/dt. A sample example with 
a = 0.25, b = 4 and 0 = 0.06 and /(O, o) = 0.2105 yields a symmetrical 
hysteresis curve. 



2.5 Ihe Statement of the Problem 


The problem can be stated now as ‘ Given the state model for the 
bar mill drive and specified rolling torque (Appendix B) obtain a 
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control configuration to achieve 

a ) an acceleration rate of O .5 p..u. to race to base speed from 
standstill with load (ref. Pig. 2-2), 
b ) an acceleration rate of O .5 p.u. to race to top speed from base 
speed with load (ref. Pig. 2 - 2 ), 

c ) a recovery time of not greater than 0,5 seconds for impact 
loading or load throw-off at top speed, 
d ) speed regulation 4 /^ 
ensuring the following, 

1 ) the armature current does not exceed 1.9 p.u. and persist 

beyond O .5 seconds duration, 

2 ) the maxim.'um rate of change of armatiire current does not exceed 

2.0 p.u./sec. and persist beyond O .5 seconds duration at that 
rate , 

3 ) the regenerative current does not exceed O .9 p.u. and persist 

beyond 0,5 seconds during the field weakening period, 

4 ) the rate of change of either of the field currents does not 

exceed 8.0 p.u./sec. and persist at that rate beyond 0.2 seconds. 

The restrictions (l - 4 ) are dictated by machine design and often 

7 

stabilising components are used in addition to the regular 
protection equipment to ensure safe operation of the mill. It 
is assumed that continuous forcing functions will be available. 




Figure 2-2 

Speed versus time characteristics. 




In arriving at the form of criterion of optimality some flexibility 

. , , . -u-i -x 15 »' 14 > 15>16 

has to he observed to ensvire physical realisability 

in optimality criterion* of the form (2.4. l) was chosen, 

(Jl 

B = min J" ji * (Njj - + B * + 1 /^ + D * (2.4 . 0 

CPhis choice was based on the following considerations J 

i ) Ihe primary variable of interest is speed and a measure of 

deviation of actual speed from desired speed should be included, 

ii ) For hard-ware realisability there shall be a specified limit on 

the value of forcing function so a function which reflects 
steep penalty as the magnitude of Y^ increased should be 
included, 

iii ) The inter-coil voltage should not exceed a specified limit 

for safe operation and this is reflected in the term 
containing Y^, 

iv ) The term B * Y^ has been included to slightly under weight 

negative values of Y^^ so as to have the peak - 10 p.u. 
reached (without raising the value of A very high). 

As such, there is no general procedure to decide on numerical 

values of A, B and B. However a convenient rule would be to follow 

T7 ' ^ 

the algorithm described by the decision table (Table 2 - 1 ). The 

condition entries are particular clock readings and maximum values of 

forcing voltage. The clocks are simulated in the program to check on I 

operational constraints described in Section 2.5.. In practice these 1 

clocks would be expected to provide protective features only. 

* It might seem that the functional E does not necessarily remain non- 
negative and convex,. As during the time Vf is negative first term c.crv\.h IvTI 
dominates in the integrand, the function E remains non-negative a 
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Table 2-1 

Key w « maxinum permissible reading any clock 
X ® reading on any of the clocks any time 
y = maximum magnitude of forcing function 
z » maximum permissible magnitude of forcing ftinction 
1 » the particular condition exists 


Condition Stub Condition entries 


1. 1 

X w !■ 1 : 1 

1 

> 

X ^ w I’ ; 

1 1 1 

y < z •; 1 ' 

> 

1 ! 

y ^ z ; 1 

L. - ■ - ‘ 

» 

; 1 

!i ; 

Action on A u increase • decrease 

'1 ! ' 

1 

decrease I decrease 

» 

^tlon on B W;; : increase i 

Action on D /' m ! ; 

j. t : 

t 

( 

increase • increase 

' 4 


The increments can be based on an expression of the form 2,4.2 to 

finally yield a value for "A” so that has a peak value of 10, 

A - A , + (10. - nJ*K(Vf ) .. .. 2.4.2 

new old ^max-old -^nax 

where K is a nultipliUr depending upon the value of 


'Viren IT-^ 
-^loV cy E 



The clocks (Pig.2-5) record the time d\irations for which the 

derivatives of field current, armature current and speed of motor 

exceeds the specified limit. The clocks are reset, to read, zero, 

only ■'nfeen O.5 seconds have elapsed after it was last incremented. 

Thus ensuring that in no one second duration conditions 

- 'A* ’B’ and 

existed heyond O.5 seconds. The numerical values 01 » 

'P' are so chosen that mill motor will not trip under normal 
operating conditions. 

The finally accepted values of A, B and D which ensured normal 
operations are given in Table 2-2. 


Table 2-2 

Finally accepted values of A, B and D. 



Value of A 

Value of B 

Value of B 


Gen, field control 


J 

1 

f 

t 


a ) zero to base 

; 4800 

: 1 

; ® 

♦ 


b) at top speed* 

1 4800 

i 0.25 

. . -t - - • ^ 

; 200 

1 

t 


Motor field control 

; 8800 

» 

: 5 

4-^ 

0 

0 

_ 


* For both impact load on and load throwoff. 



Rbutine 































CHAP!IER - III 


OPTIMISATIOg PBOCEPUBE 

Ihis chapter is devoted to a review of some known optimisation 
techniques associated with dynamical systems and suggests a new 
technique which is particularly suited to the problem stated in 
Chapter II. (equations 2.2.1 to 2.2.5). 

The new technique suggested requires replacement of product terms 
by a weighted linear combination of both the terms (involved in 
the product) with time varying coefficients. In initial guess of 
the time varying coefficient may be based on an analog computer 
solution corresponding to an admissible control. An improvement 
on this approximation is obtained by taking Gateaux differential 
of a positive semi-definite function of state variables along the 
guessed trajectory. The optimal solution is obtained by using 
parametric expansion technique in each iteration. As a result, 
one solves a sequence of optimal problems converging to the 
required solution. 

The new technique, suggested here is simple in structure and may 
be adapted for general trajectory manipulations required in 
optimisation problems. 

Some theorems concerning the new technique are discussed towards 


the end of the chapter. 
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3.1 Review of Some Optimisation Techniques 

The following is a "brief review of some of the algorithms idaich have 
"been used to solve optimisation problems. The following classification 
of optimisation methods has been observed in discussing these methods. 

a ) Gradient techniques based on iterations on boundary conditions to 
optimise a function of terminal state, 
b ) Gradient techniques for trajectory manipulations based on Taylor 
series expansion of optimising functional around the value in 
previous iteration (also called relaxation procedure), 
c ) Gradient techniques based on Taylor series expansion of state 
and adjoint variables around that of the previous iterations 
(also known as Generalised Rewton-Raphson Method), 
d ) Conjugate Gradient Techniques. 

e ) Methods based on an approximation of optimising functional as a 
polynomial in state variables with time dependent coefficients 
(also known as Paramet ric ESxpansion Technique). 

a ) Ter m inal State Optimisation Procedure^ 

Cannonical Problem Statement : 

Given x = f(x, u, t) ^ 

x(o) = Xq , t € £o, t] ] 

minimise E =« /(x(t),T) choosing u € U ^ t j *’ ** 3.t»1 

X is a n X 1 and u is a m x 1 vector \ 
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The algorithm 

The basis of the algorithm is often explained by associating some 
geometric concepts with the solution space of state equations 
It is considered that there exists a path in solution space, emanating 
from the trajectory corresponding to an initially chosen admissible 
control, which approaches that element which minimises the functional. 
This path is referred to as the path of steepest descent. The algorithm, 
therefore, reo^uires a sequence of control programs to be generated by 
calculating the tangent vector to the path of steeliest descent at each 
step. 

In general a numerical procedure will require the following computation 

1 ) Integrate the equations of motion in forward time with an admissible 

control program, (if the terminal state is fixed, one may carry a 
measure of goodness in optimising fvmctional detailing the difference 
between desired state and final state), 

2 ) Using solutions from (l) integrate adjoint equations (5. 1.2) 

back-ward in time 

p - - (v^ ^V'’\ J 

[-‘ l ) ) 5.1^ 

p{T) ) 

5 ) Determine numerically the influence functions as follows 




4 ) Modify the control program in (l ) hy using the following algorithm 


+'■;( ^) / P \ I 






i- 

) 


5.1.4 


whereCK is a suitable scalar ) 

5 ) Now repeat the steps (2 to 5 ) until x(t)) is minimised. 


Remarks 


The algorithm can be used, in general, to minimise real valued convex 
fTonctionals generated over normed linear spaces . T/ftien minimising 
functional is convex and is defined on a function space generated by 


system of ordinary differential equation, one can use this algorithm. 

The convergence of the procedure can be proved^^ when ^(T, x(t)) has at 
least two continuous Prechet derivatives for variations in u. It may be 


remarked here that convergence rate slows down as the optimum is 
apphoached. In most cases a relative minimum can be achieved. 


"10 24. 2S 

b ) The Relsixation Algorithm ' ' 

Cannonical Rcoblem Statement • 

Given x = f(x, u, t) 

x(0) * Xq , t (f 0, T 

•T ' ' 

minimise E ~ q \ ^(x, u, t) dt + x(t)) 

subject to constraint u £' U T t€ [o, T_j 
X is a n X 1 and u is a m x 1 vector. 


) 

) 

) 

) 

) 

) 

) 

) 

) 


1 ^ 1.5 


The algorithm 

The algorithm employs a Taylor series expansion of the return function E 
in. each, iteration* la general the approximation may contain higher 
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order terms , but for the description of the basic structure of the 
algorithm one needs to consider only up to the first order terms;. in. 
the series. 


The follow, ring forms the basis of the algorithm 
TTii+l /Tni+I , _i+1 ... 


= \ + P' 


(f^^^ - dt + x(T)) 5.1.^* 


In 3 • "I • 6 we , 


= / + Q. * (x^'^’' - x^) + Q F^/^u^) * (u^’*'"' - u^) I 

= f " + Q , fV ^ x ") - - x ") + 0 j 

= p^ + Q pV.-, x^) * - x^) I 


The substitution of 5-1 *7 in 5 •'1*6 results in following approximation 


• i-i j xm . • 

^ + { z: s^(t) * w^(t) dt 

0- j=1 ^ J 


5 . 1,8 


where 


=^P^(t)/C-) u. + p^. n. and ) 


1 1+1 X 

?r . = U . -- U . 

0 0 0 


3.1.9 


The variation in control signal in equation 3.1.8 is chosen so as to 
ensure the cond.ition 5.1*10 




5 . 1.10 


This leads to the following algorithm 


i+1 i .-i / i% 

^0 == ""0 " '-^0 ^ 


5.1.11 


where is a suitable constant consistant with the constraints, of u(t). 

*^0 


* Transposition is implied where ever necessaiy. 
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A typical nxjmerical procedtire would require the same steps as followed 
for the steepest descent algorithm except that influence functions are 
.obtained as described in 5* 'I* 9 . 

Remarks 

10 

Merriam discusses equivalence of the above two appro^aches when 
only first order terms are included in the Taylor aeries 
convergence rate improves considerably when second and higher order texms 
are also included. Like the first method this one also can be xised 
optimise -a real valued convex functional defined on spaces generated 
by ordinary differential equations, 

21 2 *5 2 ^ 

c ) Generalised Hewton-Raphson Rrocedure ’■ ' ' 

Cannonical Rroblem Statement : 

'Same as in equations 
The algorithm 

The basic principle involved is that of contraction mapping. So one 
attempts to obtain a transformation which ?irill map sets of functions 
into improvised ones till the mapping repeats itself. 

Define H = fCx, u, t) +<p . f > , ) 

then p-~>H/^x , p(T) x(T))/^x(t) < 

) .. 5 . 1.12 

X = ^H/^P , x(0) = Xq ) 

and ^ ~ ® along optimum trajectories ) 

From the condition <3 u = 0 one evaluates vector u as a function 
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As a first step- it is necessary to augment the set of adjoint 
equations Tsith state equations as in equation 3.1,15 

y» t) * y + gCy, t) = 0 . 

) •• 3 . 1.1 

where y^ - y^ = y^^^ » y^^ , ^ 

The algorithm requires an iteration on initial values of trajectories 
for y^.....y2^, which, to start with, may not correspond to an 
admissible control. Iterations are carried on until for further 
iterations the same boundary conditions are repeated. 

The algorithm is derived as follows : 

t) - G^(/, yb t) + (3 y) , (^1+1 _ yl) ) 

(pG/O y) * - y') - 0 

Prom 3*1 *15 it follows that , v y = 1 

On substituting this result we have the following 

* g^(y, t) + (r'gV^y^) * - y^) 

where g^ is defined in equation 5»1*15. 

Equations 3.1.15 are a set of linear differential equations #iieh 
can be solved by any of the conventional methods. 

A numerical procedure based on Newton-Raphson iteration requires 
following computation * 

1 ) Assume an arbitrary solution y(t) satisfying some boundary conditions. 

* Transnosition is imnlied where ever necenao'^.. 


I 5.1.14'" 


in-8 


2 ) Calculate the partitioned transition matrix for the homogeneous 

linearised system corresponding to 5* “I* “1 5 setting the following 


initial conditions 


y(0) - (O, 0, ... , 0, = 1,0,. 

vdiere j takes values from n+-1 to 2n. 


0 ) 


5.1.16 


5 ) Generate a particular solution y(t) equation 3 .1 •‘15 with 

initial conditions as computed in previous iteration 

{7,(0), 72(0) yjo), 7^+-,(o). ’•''•U 

where superscript i refers to previous iteration number. Let this 


solution be denoted by y • 

Jtr 


4 ) !Ehe complete solution of the non-homogeneous system 3.1.15 is 
generated by a linear combination of solutions obtained in 


steps (l) and (2) above 
2n 


i.e. y 


j=n+1 


w^ * y^(t ) •+ 


5 . 1.18 


where coefficients w. are obtained by solving n algebraic equations 
corresponding to boundary condition at terminal time. 

5 ) For the next iteration the initial conditions are chosen as 
(y-,(0), ygCo), ...» ''2* ^an 

and steps (2 to 5) are repeated until the coefficients wJs in 
5. 1.18 are small and a measure corresponding to (y ~ y ) 


approaches zero. 


Remarks 

n^his is an indirect method unlike the two methods discussed earlier 
For linear problems the convergence of the method can be tested. 


Ill 


The authors cla.im advantage in teims of computation time over 
steepest descent method because of quadratic convergence^ 

27 28 

d ) Con.iugate Gradient Technique ’ 

Cannonical Problem Statement : 

Same as in 3* *5* 

The algorithm 

The basis of the algorithm is same as that for the steepest descent 
algorithm except that by carrying the- information on gradients of t'wo 
previous steps, instead of one, convergence is accelerated. 


A typical numerical procedure will require following computations * 

1 ) Integrate the equations of motion 3 •"'•5 in forward time with 

admissible control program. 

2 ) Using solutions as obtained in step (l) integrate adjoint 

equations back-ward in time 

(p"" * ] 3.1.20* 

with p(t) = (5 x(t))/0x(t) ) 

3 ) Determine the gradient functional as 


and from 5.1.21 generate the norm 3*1.22, 

r" 

0 - 


5 . 1.21 


!'■ = - ( (g^(t) . g^(t))dt .. .. .. 5.''*22 

0^ 

4 ) The new control program is obtained by following the routine given 
in the set of equations 3.1. 25> 


* Transposition is implied tdiere ever necessary. 



« • 


5.1.25 


i+1 i i i 

U = "U + S 

where 

= *“ g^(t) + 

and 

/x^~^ = (r^ / i^ 


) 

) 

) 

) 

) 

) 


5 ) Repeat the steps from (2 to 5) "until the trajectories are 
optimised. 

Remarks 

It is a direct method like steepest descent method. It may be remarked 
that the first step will have to be generated using steepest descent 
algorithm. 

e ) Parametric Ibcpansion Technique' **^’ 

Cannonical Problem Statement : 


Same as in equations 5.1.5 except that 
E =.. 

0 



5.1.24 


Ihe algorithm 

2 

Corresponding to the set of equations 5.1.5» Bellman’s functional 
equation may be inritton as 

min ! P +^5 eA t + (Qe/Ax) . x) I * 0 .. .. 5.1.25; 

ut, R L ^ 

The method is based on the following approximation, 
n n n 

E = E(t) + SI \(t) * x^(t) + ^ (t) ^ 5.1.26 

m=1 m=1 pwl: r 

where K’s refer to some auxiliary variables. 
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If ? is quadratic expression in x and u, then it can be shown that 
one needs to consider only up to quadratic terms in expansion 5.1.26. 

This expansion was first suggested by Merriam"*^. It is easy to show 
that on substituting the series 5.1.26 in 5 . 1.25 and using the optimality 
condition one obtains the following 
u * a bilinear expression in x’s and K's ^ 

i = g(K, t) mth K:(t) = 0 ) *’ ’* 5 . 1.27 

The expression for u may be substituted in 5-1.5 so that 

X = i (xy Kj t) .. .. .... ,. 5.1.26 

The flow diagram of a numerical procedure to solve optimisation problem 

using this technique is given in Pig. 5-1. 

Remarks 

Systems governed by time varying differential equations can be optimised 

by using this method, if a finite parametric expansion for minE yields 

un-coupled initial value problems. The main advantage of this method 

is in terms of the computation time (as no iterations are involved), at 

the cost of additional storage. The forcing function u, as obtained in 

5.1.28 generally turns out to be a continuous function over the time 

interval | 0, t1. 

J 

5.2 The Proposed Tra.iectorv Ifenipulation Technique 

The optimisation procedures discussed in earlier section seem 

. fo'.. 

inadequate for the problem under consideration, because of computational 
difficulties involved which will be evident from the following example. 




Solve state equations forward in time using the structure of 
control law and stored K trajectories. 


Figure '3-1 

Steps involved in the solution of optimal -control problem hy using 
Parametric Expansion method. 











Ebcample 


Let it be required to minimise 


,T 


0 ./ 




where is specified over the time |^, 
Subject to the following constraints 

x^ = + a 2 * X2 + * X2 * J 2 ^(xj) 


* x^ + a^ * X 2 + a^ * x^ * ^(x^ ) 


^2 * a 


xj = ^ x^ + ag * X 2 + a^ * x^ * X 2 + b^Q ^ u 


u(t)CU Y t in[0, t] 


3 . 2.1 


5 . 2.2 


The corresponding set of adjoint equations are given by 


PO 


= 0 


= -2 * c^ * (x^ - x^g^) - - (a^ * ^2'^} 


Po = 


P, 


-(a2 + a^ * J2^(x^)) - a^ * Xg - (ag + * x^ ) 

-a^ * (d» j2f(xj)/(^xj) * p^ - ag * 0 jZ!/0x^) * Pg 




5 . 2.5 


The equations 5*2.2 and 5*2.5 are coupled sets of non-linear differential 
equations involving products of state and adjoint variables. Thus on 
using, steepest descent or conjugate gradient methods, one is faced with 
a formidable two-point boundary value problem. 

One may attempt to solve this two-point boundary value problem by 
solving an initial value problem by using Generalised Hewton-Haphson 
technique. The difficulty with this approach is that there is no way 
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of testing convergence of the algorithm for the equations of the 
form 5.2.2 and 5 -2. 5* Therefore, this method cannot he used until 
one has a reasonable approximation to optimal trajectories. 

For the problem at hand the parametric expansion technique cannot be 
used, because vith a finite number of terms in the Bellman's 

functional equation is not satisfied (as the powers of x do not 
balance out). Infinite series approximation is not feasible 
computationally. 

The relaxation procedure and related second variation methods of 
synthesis (based on Taylor series expansion) are valid over small 
regions of state space around optimal trajectories and hence inadequate 
for the problem under consideration. 

The difficulties mentioned above provide sufficient motivation to 
evolve a new method which should be computationally feasible. The 
proposed method requires decomposition of original optimal control 
problem in two parts. The logic of this decomposition is simple J 
one starts with a solvable control problem, which may be based on an 
approximation of system modelling equations, and follow an algorithm 
to generate a sequence of system models whose solutions converge to the 
original system’s solution. Since each problem in this sequence has to 
be solved to generate the next approximation, the original problem 
iTOuld be solved -when a specified convergence is achieved. 

The first approximation to state model is obtained on replaciaag all 
terms involving products of state variables by an approximation 
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of the form 

vhere + Wg- = 1 , 0 <i!!rj ^<'1 , 0 <rwg-<ri and 

and x^g is an approximation to x^ 

0(.^-f„) is a approximation to jZ^(x.) 

J& 3 

^ is a continuous one-to-one mapping. 

%• specifying the time trajectories x. , V i which appear in, 
the product terms, and replacing the product terms in state eqxiations 
hy the rule 5 . 2 . 4 > one obtains essentially time varying non-linear 
differential equations. As an example the equations 3 . 2,2 may be 
approximated by the set 3 . 2 . 5 > if ^ is a continuous one-to-one mapping. 

X., = a., * x^ + ag * Xg + 0.5 * * Xg * ^(x^g) + O .5 * * X 2 g * ^(x^) ] 

X 2 = a 4 * x^ + a^ ^ X 2 + 0.5 # ag * x^ * ^(x^g) + O .5 * * ^(x^) ) 

x^ = a.^ * x^ + ag .. X 2 + O .5 ^ * x., * X 2 g + O .5 ^ a^ ^ x^g . X 2 ] 

5 . 2.5 

With the modelling equations 3-2.5 replacing equations 3.2.2, the 
minimisation of the functional 3 - 2.1 can be performed by any of the 
standard optimisation procedures. However, the original problem will be 
solved only #ien it is ensured that 

I * f^3g-b| ■° •• 

Thus on employing the proposed technique the following two problems 
have to be solved in each iteration. 
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Problem 1 

Optimisation of performance functional E, mth. those equations which 
incorporate the latest approximation of the form 3.2.4. 


Itoblem 2 


minimise x, 


i - I i 


Algorithms suited to Problem 1 have been discussed in the previous 
section. An algorithm for Problem 2 is given below s 

1 ) Generate a non-negative real valued convex functional 3 •2*9 

(/(^. X ))by solving equation with the approximations of previous 

O 

iteration. 

2 ) Obtain Gateaux derivative of X ) along each of the guessed 

S 

trajectories. 

3 ) Modify the guessed trajectories by following the rule 3*2.7 

Vi .. .. 5.2.7 


k p * T) 

X • = X • — t * -i-r' V . 

•^ig ^ig G 1 


where k refers to iteration number, 

is the Gateaux gradient functional along i-th 
guessed trajectory, 

^ is a suitable constant. 

4 ) Repeat the steps 1 through 4 until 


0 T i 


and 


3 . 2.8 


A flow diagram,, detailing a numerical procedure, which uses 
parametric expansion technique for Problem 1 and the algorithm 
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given above for Problem 2 is given in Pig. 5-2. It may be noted that 
the particular convex functional is generated by an expression of the 
form 5 ♦2.9 

’^(t) = ■ ^ >_ (x. * X . -X. *x. )agr- , t[ 5.2.9 

Q_; £ j -'-g J JS 1 - J 

Some theorems and definitions, be -'ed upon vrhich the flow diagram in 
Pig.5"2 has been evolved, are discussed next. 


5.5 Some Definitions and Theorems 

Definition^ ^*-1 ( Convex function ) s Let X be a vector space, Y the 
set of real ntimbers then a mapping V defined on a convex subset 
C ^X, to the set Y is said to be convex if 


V(w * x^ + (1 - w) *X 2 ) w *■ 7 ( 5 ;^ ) + (1 - w) * vCxg) j 


7 x^ , Xg 6 C 
7 w, 0 ^ w 1 


) 

)• 

) 

) 


5 . 3.1 


Theorem - 1 ; Let S be the solution set of the ordinary differential 
equation 5.5.2. 

X = f (x,' Xg) 

where x, x C, a convex subset of S, 


3 . 5.2 


Then for each t ^ tJ 


7(t, x» X ) * ( 

g Q ) 




5 . 3.5 


1 j 

is a convex real valued, non-negative functional on C. 

Proof : For x^, Xj, C we have 

2 


(n * ^ig - * ng> 


0 


3 . 5.4 



I V — ” V.. 

j Generate first guess for state variables in product terms. 


H/ 
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therefore {.x. * x. ~ x. ^ x. ) ^ 0 

t ^ 2 

therefore \ T F' (^-i * ~ < 

o) i ^ ^ jg ^ ^■ 

therefore V is a real valued non-negative functional, 


3 . 5.5 

5 . 3.6 


To prove convexity we have to show that for any t <=- ^0, tJ and 
for any x^^ , x^^^, x^^, x^^^ , x^ 2 > ^ig 2 » ^j2» ^jg2 ^ 

v(t ) - J ^ (w * (x. ^ ^ X - x^g^ ^ x^^ ) 


Oil 


+ (1 - w) « ^2 * *jg 2 - ’^iga * ■! 

t 2 

j S » Xj^) d<7- 


0 ' i j 


I * '=^12 * == jg 2 ■ " 1«2 * * 32 ^ 


To prove inequality 3*3.7? proceed as follows : 


Let 

^ik 

**igk 

1C. *• 

igk 

""ok 

= a. ., 
ijk 


• « • • 

if 

^ii2 


then 

for 

0.^ w 



w * 

^^-ijl ■" ^ii2 


j1 ■ 

^ij2 



or 

W * 

K 31 ■ 

Hj2^ 

+ 2 

^^i02 

$ ® 

ij1 ®'i 32 

or 

2 

?r 

* (-131 

■ ^io2 


2 * w ' 

* a. 

ij2 

* (a. .. - a. ) 
ijl 132 '^ 







W * 

('2 2 >, 
K 0 I ~ "i 32 ^ 


(w * (a^^^ - ^ ■*■ ^ij2^ ^ ^ * ^±j1 " ^ * ^ij2 ^ij 

^ n 2 2 , s . 2 


5 . 5.8 


or (w * a^.,| + (1 - w) * ^ ~ * ®'ij 2 5.5-9 





^ij2 ^ ^ij1 could start with the inequality 

(1 - ») * (aj^j2 - ^ (a..2 - a.j,) .. .. 5.3. 10 

and arrive at the result 5 *3 *9. 

From 3- 3*9 it follows that for each t <r £o, tJ 

j 51 Z. (w * + (1 - w) * 2l Z1 W * (a£^.j) dQ- 

t p 

+ ( £ 2: (1 - w) * (a. .„rd(r'.. 3.5. 

0^ i i 

On substituting the values for ^^nd 3 .3 *8 we have 

from 3.3*'''^ 

,t 

J h (’'*(’'11 * ’' 1*1 - ’' 1*1 * ’'ji) 1 - (1 - ») . (nia . ” 1*2 


3.5.11 


O' 1 3 


- "igj * 


<0 ! * ^1«1 - ’'igl * ^'jl’ 


■ 2 
+ \ “• w ) * ^^i2 * ^ig2 *" ^ii2 * ^j2 ^ do 5.3.i2 


0' i 3 


The inequality 3.3.'I2 establishes convexity of the functional T for 


each t. 


Remark - 1 : V is a continuous function on C for each t. 

Remark - 2 i V has a lower bound eqxxal to zero and V = 0 -vdien 

either x., “ x. ' ^ i, for all k or x. , “ 0 Vi 

. igk ^ ^ 

Definition ~ 2 ( Gateaux differential. Gateaux differentiable, 

Gateaux derivative) ? Let X be a vector space, Y a nor mad space, 
and V a mapping defined on a domain A ^X, to the set BoY. iHien 




for any x A and arbitrary h ^ X, the limit 5-5»'13» 

Lim I V(x + a * h) - T(x)j/ a ) 

a^ 0 ^ ) .. .. 5 . 3.15 

inhere 'a* is a scalar. ) 

if it exists, is knoim, as Gateaux differential of V at x with 
increment h, and is denoted as DjjV(x, h). 

Further if the limit exists V h X, the mapping Y is said to be 
Gateeu x differentiable at x. 

If V is Gateaux differentiable then Gateaux derivative of Y at 
the point X is defined to be the mapping : X-»Y given by 

V'(h) » Dg(v(x, h)) ... .. 3 . 3 . “14 

Remark - 1 2 %en Y is the real line, then Gateaux differential of V 
at X with increment h, if it exists, is given by 

D„V(x, h) = d (V(x + a*h))/da | .. 3.3.15 

b j 

a=0 

Remark » 2 s There are three terms defined above 

1 ) Gateaux differential 

2 ) Gateaux dif.ferentiable 
5 ) Gateaux derivative 

Ehch term has to be used in appropriate context. 

Theorem - 2 : Let V be a real valued Gateaux differentiable ftmction 
on a vector space X. A necessary condition for T to have an extremum 
at Xq € X is that h) = 0. 
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I^oof ! As ^ is Gateaux differentiable and has an extremum -ai 

Xq <f X, it follows from ordinary calculus, that Y h6 X 

(d (y(x + a *h))/da) 1=0 .. .. 3.3.16 

i=o 

Remark ^ 1 '• Let S be the solution set of an ordinary differential 
equation 3 . 5 . 17 > 

i = f(x, Xg) .. .. 3.3.17 

Y a convex real valued Gateaux differentiable 
functional defined on a convex subset C ^ S then a 
necessary condition for V to have an extremum at 
Xq 6- C is that I>gY(x, h) = 0 ^ Y h e X. 

Theorem - 3 • '^or the conditions given in Theorem - 1 , the function Y 
is Gateaux differentiable. 

Proof ! To show that V is Gateaux differentiable, it is sufficient 
to show that the limit of the integrand in 3*3 exists, for each 
t 1^0, T~| , for arbitrary h Cm S, 


Let 'a' be any scalar then for any h- , h^ ^ S we consider the 
following limit : 


Lim 1. A- |i_^^^i + a.*h^) * (x. + a*h. ) ~ (x^ + a * h^) * 

a~v 0 i 0 L“ ° 


(xig ^ - (^i 




* =^36 ' ^3 * 


3 . 3. 18 


“ Z-' E. Lim 
i j a-e^0 


r ( (x. * X . + a^ * h, * h . + a * h. * X . + a * h . * x. 

L 38 i 38 i dg 38 ^ 


-x^-Xig-a *hig-h.-a*h^.x. 


- a « h 3 » xjg))‘ 
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Let X. * x.g - X. ^ x.g = h. * h h . = q^i.i 


jg ^ig" “o ^2io 


h,. * x.„ - h,. * X. = q.,. ■> h. * X. - h. * x. = q.. . 
1 Og ^ Ig 5ij Jg 1 Ig J 


Then 5* 5- "I? shall be 

■ f e'Jo ( ■"" * SiJ * 


■ I Stj ‘luj 


2 1 / 
ilij - / ^ 


E iL (h. * X. - h . * X . + h. * X. - h. * X.) * 

i ^ ^ 3S 3 ig Jg 1 ig d 


5.5.20a 


(xj * X ■ — X • -x X • ) 

1 3g 3 ig 


5.5.20b 


Expression 5.5.20b demonstrates that the limit exists Y h tr X. 



Theorem - 4 : Let S be the solution set of the ordinary differential 
equation 5.5 . "ITj ^ s- continuous, real valued, convex functional 
defined on a convex subset C'SS. Then 
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1 / i X 1 N - 

^ 1 g 2 ' *'•’ ^ 1 gn^ ^ 


5.5.21 


implies 

= Xg - ''■7(jV(xg, h) * a * h e C 

yrhore l) h = (h^, h2> ...j 
i 


5.5.22 


2) V h) * a * h = (%V(x^g^. * a * hy ... 

. • • r , h ) * a * 'h ) 


. "A ' » n 

« 1gn n 

foj isny rc‘''l 'a' satisfying the inequality 5 . 5*25 

0 aj ( V h) x- h || 5 ie) 


5.5.24 


There S(^) refers to a neighbourhood of for a corresponding £ 

§ 

neighbourhood of 'V(x^) Trith { 0 as determined by continuity of V 

g 

i 

in , 

Proof ; Let Xg f- C, Y is a continuous function of x^ defined on C, 

i+1 


For any x^ C- C and f:/0 Y^l ^ ) 6 K(v(xi), £■ ) 

& & 6 

We have = V(xi) - VV * - x^) + 0 . .. 


S 


s 


8 


where ^ Y = (<3 Y/^^ Xj^^ , 3 y/^ x^g2 ’ • • • » 2> 
and Xg'*’"' £ N(xg , ,^(i£' )) 

“^gii c & 


5.5.25 

5.5.26 

5.3-27 

5.5-28 


So now if one follows the rule 5.5*22 ensuring 5-5.29 

H *h.,, ..., * a ^ 5.5.29 


‘1g1.' 

Then Xg"*”"* £ C. 
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It follows froni 3«5«29 that th.6 scalar ® is to h© chossn 
ensuring 3.3.5O. 

0 ^ a > li h) . h || < g (t ) . . . . 5.5. JO 

where S (^ ) refers to a reighhoTirhood of x for a corresponding S 

o 

neighbourhood of with (-> 0 as determined by continuity of V. 


Theorem - 5 * Given the conditions of Theorem - 3> a necessary 

k 0 

condition for the sequence x. to converge to some x. is that for 

"38 38 ■ 

a fixed scalar 'a' and increment h 




5.5.51 


where k refers to iteration number. 


This result is an immediate conseouence of previous theorems. 


Remark - 1 s The process of convergence will slow down as iterations 
proceed, because for fixed 'a' and ’h' we notice that 
step sizes reduce as iterations proceed. 

Remark - 2 : For the V function defined in 5.5.5 the sequence will 
converge when 



as this is the condition Tidien = G. 


5.5.52 
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CEiiHm ~ IV 

THE COITJfiOL 0PER ^ .TI0NS QU MlilS DRITO 

In the previous tvo chapters the modelling, and optimisation procedures 
were discussed. In this chapter the problem stated in Chapter II iitll 
be solved by using the optimisation techniques discussed in Chapter III. 

It may be recalled that during the generator field control mode of 
operation the state eouations are linear and during the motor field 
control mode of operation, the state equations are non-linear and 
involve product and functional types of non-linearities. For both 
the sets of equations parametric expansion technique has been used to 
obtain tho control program. The use of iterative trajectory manipula- 
tion technique described in Chapter III is demonstrated in computation of 
a control law for motor field control mode of operation. 

This chapter includes discussion of computational aspects like step 
size selection and selection of integration procedure etc. The 
resTolts of computation are given in various graphs and tables. 

The chapter ends with a description of a hypothetical situation 
detailing computer aided execution of optimal control program for 
operation of all speed ranges of the main drive. 

4,1 Speed Control by Generator Field Voltage 

The problem of speed control during generator field control mode of 
operation -can be stated as follows s 
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Given the following state equations 

d Ij/it . + 1,1 * Tjg 

d I^dt - aj, ♦ Ijg + e22 * + =>5^ * \ * 

d V« ■ ^42 * ^>2 * ’’l 

d Ifjj/dt = 0 , = 1.0 p.u. 

where ® ~l3-( ” -0,288 

^■21 “ -®22 “ "^24 “ 

= h^ — 2 .24 


4.1.1 


Find y^g to minimise 

T “7 

E. J. [ A . (N^ - + B . Vjg . . B . (V^ - dt 4.1.2 


where and refer to actual and desired terminal voltages, 

V, - K * I. - r * I 
■t g fg ag a 

\d ■= ^ when rolling takes place at top speed 
D * 0 when reversals of speed take place 


Ihe magnitudes of A, B and I> are obtained by following the logic 
of iteration discussed in Chapter II, so that, the following conditions 
are ensured t 

1 ) The maximum magnitude of does not exceed lO.Op.u. 

of is continuous. 

5) jRie performance specifications discussed in Chapter II are 


2 ) The form 


achieved. 
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The Hamilton— Jacobi equations corresponding to this problem can be 
written as 

\2 . ^ . „ . , .. s 2 


f“.. ^ ^ + B ‘ Tjg + B » (Vt - f 3 E/2)‘ 

Vfgfc ^ 

C,Wa Ifg) * + 0^d\) * \ j 4.1.3 

Let E = K^(t) + K^Ct) * + K^(t) * + K^(t) * + jEf^(t) * 

+ KgCt) * Ij.g . + K,(t) * KgCt) . 4 + Kj(t) * . N 

*• ♦« ■• 


*hoM 


On substituting the expression 4* "I -4 and equations 4*1 •'' in 4- "I -3 
one obtains the control law 4* I* 5 and auxiliary differential equations 
4.1.6 as follows; 


Vfg . - b, MKg 2 - + Kg - la + - Nj - B 


4.1.5 


K/dt 

+ 

^4 

* b2 ■> 

.Tl ^ 

I- A * - 

md 

0.5 

(K 2 

* b., 

C\J 

+ 

= 

0 

i^dt 

+ 

s 

* ®21 


* b2 * Tl ■ 


"S1 

- 2 

* K * 
5 

b. 

* 

I 





(K 2 1 

b.j + b) = 

0 






Ky/dt 

+ 

"S 

* &22 


* a^2 * 

+ K, 

* bg 








-*'6 

‘ b, * (Kj 


+ B) 

= 0 




K/dt 

+ 

s 

* ^24 

* + 2 * 

So * 

• bg * 


2 * A 

* 

\d 






* ^1 * 

* b., 

+ : 

B) = 

0 



ly'dt 

+ 

2 ^ 

^ * 

^11 

1- Kg * a2., • 

- 2 

k2 

.b^ 

= 0 



Kg/dt 

•f 

^6 

* ^11 


* a22 + 2 


* &2^ 


* ®'42 

* 



2 * * Kg * b., * 0 
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d ydt + ^6 * ^24 * + ^9 * ^21 ) 

- 2 ^ *“>1 = ° I 

d Kg/dt + 2 * Kg ^ a 24 * + K^ * a ^2 ^ 

-0.5*K^*-5^ =0 

d lydt + 2 * Kg ag^ * ^9 * ^22 ^ ^ * K^q * 

^42 * - ^6 ^ ^ * ^1 "0 

d K^o/dt + K^ * *^(I^) + A - 0.5 * 1^ * Td^ =0 ) 

k^(t) - k^(t) «... - k^q(t) -0 

To obtain the optimal solution one first solves the auxiliary 
differential equations 4 * "I *6 back-wards in time and stores those 
trajectories •"Jiich are used in obtaining optimal control law. Now 
using the control law 4*’^*5 one solves the state equations 4 . 1.1 in 
forward time to yield optimal state trajectories. The discussion on 
computer results is given in Section 4-5 after evolving the computational 
strategy for the motor field control mode of speed control. 



4.2 Speed Control by Motor Field Voltage 

The problem of speed control during motor field control mode of operation 
can be stated as follows » 


* N 


m 




Given the following state equations 
d Ifg/'it - 0 

dl/dt - agi +a^2 * \ ^ 


) 

) 

) 



d Ij/dt . Ojj ♦ » Ifi. + b, » 

a ft/dt . a 2 . ^(I ) * 1. + b * Ti, 




"fm "ben ,<0.75 P.u. 


0.75 + 0.8 * (I^ - 0.75) 

0.95 + 0.1515 * (Ifjjj - 1.0) 


4.2.1 


•• •• (adjusts field time constant) 


= -bj = -0.218 


Find. V to minimise 
fm 


E - ( [ A . . B . A + B . (V^ - ] dt 4-2.2 


where and refer to actual and desired terminal voltages. 


) * N + I * r 
^ ^ fm m a am 


V+,:, =1.0 p.u. 


The magnitudes of A, B and are obtained by following the logic of 
iteration discussed in Chapter II, so that the following conditions 
are ensured. 

1 ) The maximum magnitude of does not exceed 10 p.u. 

2 ) The form of is continuous. 

5 ) The desired performance specifications aa stated in Chapter II 
are achieved. 


The Hamilton-Jacobi equations corresponding to the above problem 
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can be written as 
min 

V . £ V 
fm 




* QV,'. TJ . H- 01;^ I^) . \) * *^[1 * “ 

4.2.5 

To use the iterative trajectory manipulation technique, the system 
equations are first approximated to the form 

4 if/at - 0 

i l/at - ^21 » Ifg *\+ 0.5 » .{N^ H^) 

4 V" - "55 * ■'fm * Ifm ^ *>5 * ""fla * 

d N /dt - 0.5 * a;, * I~ „ * I„ + 0.5 * a' * I. * I + b^ * T, 
nr '’ 4.2 fmg a 42 fm ag 2 ■ L 

v;here a ^2 modified coefficients corresponding to a 2 ^, a ^2 

to have piece-’’i 80 linear approximation to versus 

characteristic. 


) 

) 

) 

) 

) 

) 

) 

4.2.4 


^fmg’ ^ag \g ®oi'^®spond to the guesses of the solution 
and respectively. As a first guess and are 

generated from the analog computer solutions discussed in Appendix A. 

Let 3 . K, 4 4 Kj . 4 + 

■■ 4.2.5 

On substituting equations 4.2.4 ^^4 4.2.5 irr 4*2.5 and applying 
optimality conditions one obtains the control law and auxiliary 
differential equations as follows 

fe “ “ ^fm * ^5 "" S * \ * ^fm + 5/^5) 4.2.6 


V. 
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d K/dt + Kg * a2-, * Ifg + \ °*5 * 

(K^ * -b^ * K^ + B)^ + 2.5 *D =0 
d Kg/dt + Kg ^ agg + 2 * ^ + K;^ ^ bg * Tl - Kg * b3 * K^^# 


^ + b) 

+ (0.5 

*K^- 

D) 

a ^ 
^42 

* ] 

'■fmg 

a ’ * ir + K, 

24 mg 4 

* a!„ 
42 

* ; 

[ ) + 
ag' 


* a. 



B) * b 

5 * 



* b. 


* (I 

* a!„ 

+ 2 

* u„ 

*11 ) 

fg 

ag 

42 


fm 


mg 


fm 


d K^/dt + ag^ * K^ * I^g + 0.5 * Kg * a^^ ^ + 2 * (- A * 


*^2 -TL-n^-D*I^g + 0.5 

(Kj *bj +B)) = 0 

d iC^/dt +2 * K^ * *^22 ~ * ^^^6 * ^5 * ^fln^ “ ^ * ^42 * ^ftag^ 

+ ^ * (=42 * ° 


d K^dt + K^ h * ^^42 ®'35 * * ^4 


42 



(K^ * : 

^ag 

^9 " ^ftng^ - 2 

* Eg * Kg * (b^ 

* V - o 


d Ky/dt 

+ IC^ * 

^■42 

- Kg * K^ ^ (b3 

* ■" t’S * 

“24 ■ 



+ K.,q 

* ^42 

) * I =0 

^ fng 




d Kg/dt 

+ 0.5 ^ 


*^2i 

**42 **ag> * 

D * ((a^g * 



+ (''fB 

* N f)+ 2 * K„ * 

mg' 0 

■'to * Sj • ^ 

* 

= 0 
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d K^at + 0.5 * * (K^ * 

- 2 * Kg * * (b^ 

d K,(/dt + A + 14, * (0.5 * * I3 + B * 4 ^) ] 

+ 0.5 * IC * a' * I = 0 ) 

^ 1 24 fmg , 

K,(T) - K^Cf) K^q(t) - 0 

where is used as a coefficient to compute terminal voltage 

(a function of in terms of 


N 


+ * I 

6 fmg 


) + K^q * a;^ * I 


%2 


*^9' ^33 * ^rm 


= 0 


) 

a€ j 
) 


To obtain the optimal solution to the original problem (corresponding 
to equations 4»2«'1 and 4*2.2) one first solves the optimal control 
problem corresponding to equations 4.2.4 and improves upon these 
solutions till the solutions corresponding to the original problem 
are obtained. 


The atuciliary differential equations 4 * 2.7 (corresponding to the 
approximation 4.2.4) are solved back-wards in time and. those 
trajectories tidaich are used in obtaining optimal control law are 
stored. This storage is in addition to the storage required for the 
approximating trajectories ^ag^* control 

law 4.2.5, one solves the state equations 4*2.4 in foinurard time and 
also generates a convex real valued non-negative functional 
corresponding to the expression 4*2.8 (notice the similarity with 
the expression 5*5*5) 




# N .* 

mg 




. 2.8 
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To improve upon this approxinicLtion it is required to evaluate Gateaux 

Differential along the trajectories I , and IT . 

img mg 

To evaluate Gatoaux differential, one first obtains a solution of the 
set 4.2.1 for some continuous form of ^fm (preferably with oscillations), 
and use these solutions to generate variations on each of the trajectories 
and For this computation a sinusoidal perturbation 

over-riding an a.draisslble control was used to generate variations to 
be used in evaluating Gateaux differential. The frequency of the 
sinusoidal perturbation was varied as the iterations proceeded. The 
numerical procedure detailing the evaluation of Gateaux differential 
is discussed later in this chapter. 


The trajpotories are modified next by following the rule 4.2.9 • 

N - ) 

mg new % old ^ tng old op / 

^ag new " ^ag old " ^^^^ag old’ ^ap^ j .. 4-2.9 

^mg new * ^fm old "" ^ old’ ^fmp^ ) 

Tdieru tho subscript p refers to the solutions obtained with a 
perturbation riding over an admissible control to determine the 
general direction of improvement fer the guesses ^ag ^mg* 


This modification schema is carried on till 


i N - s ! + li - I I + ~ I- j < 0.0001 •4.2.10 

\ m V^mgl . j -^a *ag | ) fm fmgj 

which ensures that the original problem is solved optimally if 
^m’ ^firg/ ^agV^mg* ^fii®’ ^ap’ ^mp^ “ 
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■1 . 5 Computational Considerations 
Soleotion of step size t 

The optiEiisp'tion algorithm employs parametric expansion technique. So 
solutions to the a.uxiliary time variahles are first obtained back-wards 
in time and these trajectories are stored to be used in obtaining 
forward time solution of the state variables later. It is obvious 
that if the initial step sizes (corresponding to the terminal time 
in (O, T))aro small, than error propagation toivards the end (corres- 
ponding to initial time in (O, t)) would he small. Therefore the 
computationol grid near terminal time must he very fine. 

The analog computer solutions (feasible solutions discussed in 
Appendix A) indic>"'te hi^ rote of change in state variable trajectories 
during initial portions in time interval (o, t). This would suggest 
that step sizes near the initial time must be small to obtain the 
reauired numerical accuracy in the solutions of state variables; 

(Por all computations here, the required numerical accuracy was 
placed at sight significant digits.) 

The arguments above suggest simaller interval size at both ttie endji of 
every time slice (O, T). For this computation eveiy time slice was 
divided in 5OO intervals such that the step size decreased linearly 
towards the ends as described by the expression 4»3»1- 

hjj " ^n-1 0.00001450* ((250.1 - n) / ( 250.1 - n)) * T .. 4-3.1 

■where n refers to the number of interval and h^ =■ T / 5OOO. 
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On^ can, in general, generate an error botind^^ controlled step size, 
hO'.^ovcr, for operational considerations a fixed routine, like 4.5.1, 
is to -be preferred. Besides this, a strong justification for using 
the expression 4.3*1 lies in tho fact that on halving the step sizes no 
particular improvement in eighth significant digit is observed. 

Saleotlon of method s 

Bie variable step size restricts the choice of solution methods to 

19 

Rungo-iiutta type routines. Sill *s modified version of fourth -order 
Runge-Kutta method was used in tho present computation, Ihe round off 
errors accumulatod at each step were found to be smaller than ei^ts 
significant digit. Bie smallness of round off errors at any step 
can bo ascertained by following the procedure suggested by Gill. The 
procedure requires a comparison of the results Tiihen the intermediate 
variabl 'S are stored in double precision with those obtained idien 
these variables are stored in single precision. If the two results 
compare well than the smallness of round off error in every step and 
overall numerical stability both are ensured. 

4*4 Computer Program 

®ie complete flow diagram is divided in four parts so that different 
aspf^^ots of the program can he emphasized evenly. 

•i 

Fig. 4-1 s The flow diagram in Fig .4^1 essentially depicts the rolling 
schedule. The switches P, PF, FFP, PFFF, GG, GGG, GGGG and BP combine 
to yield different modes of operation. 



Switch P* • +1 for generator field control mode of operation. 

-1 for motor field control mode of operation. 

Switch FF J +1 for rolling at top speed of operation. 

~1 for rolling helow and at "base speed of operation. 

Switch FFP s +1 for - base speed to + base speed of operation. 

~1 for + base speed to - base speed of operation. 

K 

N^d(K) - - ^ (2 - (^ h.) / T) .. 4.4.1 

Switch PFPF I +1 for + top speed of operation. 

-1 for - top speed of operation. 

N^d(K) * *2 .. 4.4.2 

Switch BF I +1 for loading at top speed of operation. 

-1 for load throw-off at top speed of operation. 

Switch GO : +1 for field weakening condition (to race to top speed 
from base speed). 

-1 for field strengthening condition (for changes in 
speed from top speed to base speed). 

Switch OGO ! +1 for + base speed to + top speed of operation. 

-1 for - base speed to - top speed of operation. 

K 

N .(K) - GGG * (l + ( ^ hj / T) .. 4.4.3 

incL 

Switch GCKJG s +1 for + top speed to + base speed of operation. 

-1 for - top speed to - base speed of operation. 

■ E 

n Ak) ^ GGGG (2 - ( 5_. h.) / T) .. 4.4.4 

ilia • ^ 

In opeifational stage these switch combinations will generate the 
schednling algoritiim for desired speed of operation. 


Ehe tiro iwitch positions are denoted as *^1 and -1# 
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4~?. * fl<'w diagram in this figure refers to common program 

'Bhich is used during both the modes of operation. The flow diagram, 
as it io drawn, assimes switch position P as -1 i.e., the motor 
field control mode, ’'Ihen the switch position is +1 the function 
generation routine is skipped. 

The phrase "for all J" refers to a DO loop in POETKAN EROGRM, 
Similarly h^ refers to n-th interval. 

Various coefficients used in generating intermediate steps in 
Gill's routine in fact can be generated once and stored to be used in 
computing intermediate steps in solution of the auxiliary equations 
as wall as the state equations. 

During operational stage one need not store the matrix K(j, N). One 
will only store those trajectories idxich are used in generating 
control law as this is the only useful information during a particular 
time slice. 

Pig. 4-3 * The flow diagram in this figure refers to the common 
program which is used during both the modes of operation. The flow 
diagram is drawn with the assumption that switch P is in -I position, 
^^enever switch P is in +1 position the function generation routine in 
the flow diagram is skipped. 

A particular feature of this flow diagram is the simulation of 
protection clocks. The flow diagram for iiie clock incrementing 
routine is given in Pig. 2-3. Bere the clock routine is simulated 
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Perturbation : 

Tbe perturbs tions sre gensmtod. from the solution of equstion 4»4*5 


.k+1 „/ k+1 ,k+i X s 

X » fU , u , t) ) 

■"here u^”*""^ »= u^ + /J * (k ) ) 

' I ) 

/! » Sin f . 1 * k * t ) 

" ) 

,/4 “ 0.005 ) 


4.4-5 


This perturbation is used as the increment ‘a *h' as described in 
Section 5 . 5 . 


OateavDc diffc'rontial t 

To gen«.rnte the Gateaux differential for each n (n takes the values 

from 1 to 500) on each of the trajectories of I„ , I and N the 

fmg’ ag mg 

follo’^lng steps are followed 2 


1 ) Store the trajectojy (say in auxiliary storage. 

2 ) For n = 1 through 5OO, generate 

W new * ^fmg old * (corresponding trajectory in 4 - 4 . 5 ) 

4 - 4-6 

3 ) Oorreaponding to the change 4.4.6 compute ^fm new» \ new 

la and simultaneously generate ^^^(t) 


v+i(t) - , fd 


* I. 


- 1. 


* I 




Q j 1 } ^ new fmg new fm new ag old + 

f 


(I ♦K - 

^ ag old m new mg old 


a 


N * I / 

mg new new mg old fm new Jd <7- 4.4.7 


by smBming the corresponding stair case approximation of terns 

in 4 * 4 ''* T* ■ ■ ■ 
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4 ) For n * 1 throxx^ 5OO generate 4.4.8. 

^fmg new * ^fmg old " * (corresponding solution in 44.5) 

. 44.8 


5 ) Oorrcsponding the change 4.4.8 compute , N and I 

fm new' m new a new 

end simultaneously generate 4.4.9 
t 

V -(t) - I f(l * 1 ^ - I * I f + 

0 * 


a new fmg new fm new a old 

* N - N * I + 

fmg new m new mg old fm new'^ 


^ V old * \ now * V old * new^i^ ^ '*• 4.9 

6 ) 'flnie Gateaux diffcrontial is now computed as 


old’ ^fm perturbed^ ** ^^+1 " / 0.01 .. 4.4. 10 

7 ) The w.ldo of %V(I^ is storod ao>d stops 1 

to 7 aro followed for I and N 

ag mg* 


8 } 05 x 0 trajeetories are manipulated next by the rule 

4 Perturbed) / * -1 * >=) | 

i perturbed^ / (1 + .1 * i) 

N^1 . - 0.0001 . %V(»^, p^rturbed^ / (1 t .1 * b) ) 

9 ) ®ie steps 2 to 9 repeated till the inequality is 


satisfied. 
For each n 


fmg 


- I- } 
fm’ 


+ 11 - 

I ag o 


I * 1 % - "ml 


4.4.12 


* The multiplier 0.0001 limits the largest tariation in tra^toiy to 
about 5/0. The factor (l +4 * k) »akes the step size a function of 

iteration mimber. 



The inenuality 4 *4 >*12 corresponds to the condition that Gateaux 
differential of the function is miniinuin and function V has achieved 
the lov'er botmd. 

4.5 Computer Results 

The Table 4.I sximmerises all the results detailed in the Figures 4.5 
through 4-7» 

The main observation from theser esults is that hi^er loading of the 
main drive would be possible with the given specifications. 

The Figures 4*8 through 4 •18 demonstrate how the trajectories get 
manipulated inhen the algorithm described in Chapter III is used. 

Figure 4 •8a • shows the manipulation on the solution of the; speed 
trajectory. 

Figure 4»8b t shows the manipulation on the guess of the solution 
of the speed trajectory. 

Figure 4.9a » shows the'~manipulation on the solution of the motor 
field curr(?nt. 

■^igure 4.9b s shows the manipulation on the guess of tire solution of 
motor field current. 

Figure 4.10a : shows the manipulation on the solution of armature 
current. , 

Figure 4.10b : shows the manipulation on tiie guess of the solution 
of armature current. 



Try of the computer results 
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GILL'S ROUTINE 


M 



I 


For All J 

KI(J) ^ - _ h K'(J, n - 1 ) 

K(J, n).... K(j” n - 1 ) + (1 ^^" 075 ) * (Kl(J) - KQ(j)) 

K:i(j) (2 +- 42 )- *KI(J) - (2 + 5 *4 0*5) * 


; KI(J) 

I K(J, n) -- K{j,n) + Kl(j) / 6 - KQ(J) / 5 


Fourth Step 
For All J 


No 


S/ 


Yes 


n 0 


Store the 
matrix K(J, n) 


To solve for state 
equation, as in Fig .4-5 


Figure 4-2 

Solution of Auxiliary time varying differential equations. 




Figure 4-5 

Solution of State ©■'■ufltions. 




4' JilHS'j? Ai'j:'UOXlm5tlOK 




l.i".'.l c. n'-rol pr.'. lem with latest approximation'' 


i Obtain HFF - X(L,N) - X (L,F) FOR ALL L, FOR AT.T, kt, 

• r , ■ fi , „ 


Yes 

<^0.0001 



The problem 
is solved. 


Oen-'rate rond store the solution iwith a perturbed trajectory of an 
adinisaible cmtrolf The perturbation frequency is varied linearly 
an iterations proceed. 

m . ' V ' .I «m' - nM 


L + 1 


' . x(l, II) , in an aux iliary storage 


FOR ALL N. 


i I'erturb the tra jectory X (L,N) by the rule FOR ALL U 

; ^ H- e- . 0 (L,IT) 




Solve the optimal control problem for this approximation. 


Obtain 5‘. (X^(L,N) * X(J,N) ~ X^(J,N) * X{L,N)) 


L J 


L 2 J 


S 


FOR ALL N 
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‘ crturb the trajectory X (L,N) by the rule 


FOR ALL N 




~3olve the optimal control problem for this approximation and 
generate a as in step 4* 


! Ooinj)Ute Gateaux differential for L-th trajectory by the rule 
i . (V_^^ - / (2 . c- ) FOR ALL H 


Restore the perturbed trajectory 


FOR ALL N 


Ro/V Yes 
L > 



I Jfodify"x '(L,N) ^X (L,N) + D(jY(l,N) 

I S O . . I- ... 


FOR ALL L 


Key ; X(l,N)=Speed of motor 

X(2,R)=Armature current 


Figure 4-4 , . 

F1.0W diagram for tr?i.jectory manipulation. X(5,N)=Mbtor field current 
*CorreG ponding to analog computer solution, X^(L,N)sApproximation to X(L,ir) 
T^ith sinusoidal perturbation. 





Igure 4 • 50 - 

^peed s zero to base 

Operation mode? Generator field control. 







mataire terminal response 

de of operation: Generator field control mode 







^ 0/4 0/8 1.2 1.6 

time in seconds 

Figure 4 *66 

Armature and field response at top speed of operation. 
Made of control - Generator field voltage 







Pigxire 4 '713 
Field response 

Mode of operation - Motor 'field voltage control. 




Figure 4 *70 

Armature terminal response 

Ifode of operation - Motor field volfage control. 



Iterations on the speed trajectory. 




Iterations on the appa;'oximation to the speed trajectory 











Iteration on the solution of armature current. 
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\J' ' ion of 

•^ic •■ jr t'r.n ^nr 1 ic renHona w^i,y computer aided operations are 

r>'Cnr;"x j- 1 '’.e'l In practice * 

1 ) rh* execution of control lo^^ic requires initial storage of 

ar.-.rOfrinte initial conditions of the auxiliary variables. During 
the inotcr field voltage control mode, storage of the guesses of 
tl.f Htate '.n»ri'>hlo trajectories will also be required. 

2 ) Tiie tercui* wraua time characteristics will vary as the 

rnllinK prnco' in and thofefore during the handling time*^^ on 
the work piece the tornue versus time characteristics will have 
t.i. h< •,nd predicted for the next pass. Also during 

thin tim'' %h*' initlnl conditions of the auxiliary time variables 
•ill have to b'.- computed, 

tt "i?iy bi- ''( marked here that the modes of operation of the mill main 
drive control p.ystem hns follo”»ing fixed patterns 
a ) durinf! the top fl|x'ed * Generator field voltage control 
b } ii' iwHon the toj and base speed « Motor field voltage control 
c ) hetwjen b'l.'i!,’ apt’od ■ Generator field voltage control. 

Tnuij* iTiu follow with a certain regularity because of the 

nature of the rolliiv operations involved. The regularity of the 
pittextiB inakfta thti use of computer aided operation attractive 

*■ Tirii. required to {jooition the work piece or some auxiliary 

conti-olU/r. 



1, o'- tho K,oiiibUity of introducing ‘learning capabilities' 

o rono..‘^u-nc'^ of dnt^' logging facilities available with modem 

Q 

pr cona control ccin'utera like ft-odac 580 or Tosbac- 7000 . This 
ic cvpl.-in-rl in d<. tails in the Figures 4-II through 4-I5. In 
Figuroa /t-n through 4-15 it is assumed that a single central 
processor -■'11 bu r^'quired to achieve execution 0 ' number of 
progrona (siirultont ounly residing in memory) in such a way that 
none of iho progrpmo need bo completed before another is started 
or continui d. Such a facility is generally termed as multi program- 
ming focJllty, 11110 la desireable because though operations are 
8c^i‘> .ulod by tb*u m^tln program or an executive program, the computer 

should be ruf'dy i-o receive and transmit relevant Information within 

55 

acc>;ptn,blt. renponae time’". 

Figur* 4-15 det'.ila the use of data logging in introducing the 
1 jprnlng cnprlllity In the overnll system. It also explains that 
In gen'-rol , tht iynamlo program modification, if needed, ean be 

onrried out. 

In brief the corainitt-r ou#t to be a real-time'*^ system. 


* Rer'I-timo (idorliy) is tin extreme case of on lineness w^ere the 
re«}X)no^ tiiM of the computer systan shall be zero. 



h> cauue of the possibility of introducing 'learning capabilities ' 

as a consequence of data logging facilities available with modem 

ft 

process control computers like ?rodec 580 or Tosbac-yOOO, This 
is explained in details in the Figures 4 -II through 4 -I 3 . In 
Figures 4“'1'' through 4“'13 it is assiimed that a single central 
processor will be required to achieve execution of number of 
programs (simultaneously residing in memory) in such a way that 
none of the programs need be completed before another is started 
or continued. Such a facility is generally termed as multi program- 
ming facility. This is desireable because though operations are 
scheduled by the main program or an executive program, the computer 

should be ready to receive and transmit relevent information within 
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acceptable response time 

Figure 4~‘15 details the use of data logging in inti'oducing the 
learning capability in the overall system. It also explains that 
in general, the dynamic program modification, if needed, ean be 
carried out. 

In brief the computer ought to be a real-time* system. 


* Real-time (ideally) is an extreme case of on lineness where the 
response time of the computer system shall be zero. 




General arrangement of computer aided bar mi: 





t?d hJ 


Set 7fork 
Sv/itch On 




Transmit buffer of the channel indicated by channel 
pointer with the current pattern. For the channel 
pointer and current pattern combination execute the 
control policy. Compute the next channel pointer 
location. 



Set 

Interrupt 
Switch Off. 

The next work 
piece has 
arrived. The 
processed work 
piece has left. 


Yes^ \No Continue with 

All present processing 


*! 

The next work 





piece is yet 

No ^ 

-^Int errupt^ 

Yes 


to arrive. 







Figure L[-^2 

Flow diagram of the operation of the computer system. 
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Yes 

Recalculate 


./Yes 


Updating of Ifefhenia.tical model. 

Fi.gure 4-1 5s Data logging and program modification. 
'Srror Corrections in Successive Passes 


Physical Processes Measurements 
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CHAH’ER - V 


CONCLUSIONS MB DIRECTIONS FOR FDRTHER RESEARCH 

In this chapter it is proposed to cite the min restilts in the thesis 
and list some problems uihich may prove interesting for further 
investigation. 

5.1 Main Conclusions 

Ihe main conclusion from this investigation is that, with the 
experience accumulated on nature of operation of rolling mill main 
drives, it should, in general, he possible to use process control 
computer to schedule and optimally regulate the speed at ^ich the 
rolling is expected to take place. Since the control logic is 
based on a performance index tdaich has associated weightages for 
the speed, terminal voltage (essentially as a function of state 
variables and machine parameters), the protection system suggested 
by clock records of the duration of abnormal condition, is not 
only feasible but is reliable also. An observation of computer 
results reveal following J 

1 ) Speed over-shoots are almost completely eliminated. 

2 ) There is no regenerative current during field weakening. 

5 ) Variations in terminal voltage are larger than those usually 
observed during such operations. 

One does not have to get greatly concerned about the terminal 
voltage variations as regulation of termini voltage is not an 



important consideration here. This is so because of isolation of 
armatUTG terminals from the mains. This ,of-course, will be a major 
problem when one machine complex is used with heavy rectifying 
equipment connected to mains. 

This thesis is a significant contribution to the design approadi 
and technique used to solve optimal control problem associated with 
systems whose state equations contain single-valued ftmctional non- 
linearity and/or product type of non-linearity. The trajectory 
manipulation technique suggested in the thesis is general enou^ to 
be adapted for trajectory optimisation problems also. 

5.2 Directions for further Research 

The following is the list of some interesting extensions of the 
problem discussed in the thesis : 

1 ) Feasibility studies to examine compatibility of present day 

process control computers for use in speed regulation as 
suggested in the thesis . 

2 ) It will be worth -vhile to examine how the performance of the mill 

mn.TTi drive compares (with the present results), with the condition 
when one is permitted to have discontinuous forcing functions. 

5 ) An immediate extension of this thesis would be to examine the 

class of problems which can be solved by employing "the trajectory 
maniijulation as suggested. A computational comparison of liiis 
method with those available in literature will also be revealing 


in many respects. 
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APfSNDIX - A 


An iijialog Computer Study of Haskett's Sche me 

!Che main purpose of this appendix is to demonstrate existance of a 
feasible solution for the problem stated in Chapter II. The 
appendix describes a particular investigation iidiich was carried out 
to find compensation parameters for a speed regulating system 
identical to Haskett^s^. 

A.1 Ilrive Equipment and Performance Requirements 

Same as described in Chapter II. 

The schematic arrangement is detailed in Fig.A-1. 

A. 2 Tlie Logic of Regulatim System Itesign 

Though this system has been designed primarily to regulate the speed# 
limitations on machine design call for additional control loops to 
restrict the maximum arnature current and observe tolerances on 
aruBture voltage. Such systems fall in the catagoiy of multi- 
variable control systems . The logic governing the design of such 
multi -loop system is followings 

The main or outer-most feed back loop is to regulate primary 
variable of interest idiile the inner control loops are used to 
ovenfide the error signals of enveloping loop and provide a 
reference signal to the immediately enclosed loop. Thtas reference 
signal for any inner loop is obtained from the error in the 




A-2 


enveloping loop. In the system under consideration, helo-w base speed 
of motor, the error in speed serves as reference to the generator field 
current. Above the base speed, the reference for motor field ctJrrent 
is modified by over-voltage at armature terminals. The tendency of 
over-voltage Is seen during field weakening due to steep rise in 
speed of motor tshich is due to release of stored magnetic ehergy. 

A. 3 Protective Features and Stabilising Components 

The following machine protection is available 

1 ) Protective cut-out in armature circuit provided at 2,3 p.u. 

2 ) Over-voltage relay set at 1.2 p.u. 

The stabilising components are provided to achieve safe operating 
limits on, 

1 ) the maximum rate of change of field current so that punch 

through does not occur in the thyristors feeding the field 
system, 

2 ) the maximtim rate of change of armature ciarrents to prevent 

arcing across brush contacts, 

5 !) the maximum acceleration of motor to prevent breakdown in 
mechanical transmission. 

I|y limiting the outputs of both the field error amplifier and the 
armature error amplifier to preset values , maximum rates of 
changes on field ciorrent and armature current can be achieved. The 
four saturating amplifiers are shown bearing block nos. b1, b2, b5 
and b4 in Pig. A-2. Note that to obtain fast response large 
amplification is associated with these blocks. 




Figure A-1 




A. 4 M athematical Model of the System 


The per unit values for normalisation are same as chosen in Chapter II. 
With the assumptions as stated below one arrives at the block diagram 
as given in Eig,A-2. 


Assumptions, 

1 ) The generator field characteristics is linear. This assumption 

does not introduce any appreciable error. 

2 ) The generator speed is constant. Generally the ward-1 eonard 

generator is driven by a synchronous motor with over excited field, 
so this assumption is close to reality. 

3 ) The time constants of triggering components are zero. The time 

delays introduced by triggering components are much smaller than 
the effective field time constants '(effective time constant with 
forcing field is smaller than ishen the field is normally excited). 
This assumption simplifies the model. 


A. 5 Analog Simulation 


The armature circuit and mechanical circuit are simulated by 
equations A. 5,1 


(a i/dt) - / \ + (Sg - e^) / \ 


(d Ng/dt) - 


a 




A. 5.1 


The simulation diagram is given in Fig .A. 3a, b. 



The error amplifiers are simulated by a first order approximation.- 
For an amplifier Tdth gain of 5O and response time of 80 milli -seconds 
the approximated transfer function shall be 5O / (1 + p0.02). The 
corresponding simulation is shorai in Fig.A-jc. 

The simulation of Ramp is shown in Fig.A-Jd. Here we have a 
saturating block for ^diich the output is a input polarity actuated 
step. The step size can be changed by moving the arms of potentiometers 
biasing the diodes. The output of saturating block shall be of constant 
amplitude until the instant -when integrator output attains reference 
ampli-tude. As a result, the output of integrator will be a Ramp before 
a matching amplitude is reached and a constant thereafter. 

The motor field flux saturation is simulated by following a procedure 

51 

suggested in ^/hite and Woodson . The inverse saturation factor was 
calculated and a piecewise linear approximation required slopes of 
1, 1 , 3^5 and 1.657 between points 0-0. 76 p.u., O.76-O.95 p.u. and 
0.95-1 on ordinates. The analog computer simulation is shown in 
Fig.A-3 e,f, 

The complete simulation is shown in Fig.A-5g. 

Simulation Results 

As mentioned earlier there are four basic control loops, one each 
for the generator field current, armature current, speed and motor 
field current. For the first three the error in enveloping loop 
provides the reference. The compensating parameters can be computed 
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Figures A-3a,ls 

Simulation of armature circuit (A-53')» Simulation of Mechanical Circuit (^- 5 ^ 5 ). 
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Figure A-5 c 

Simulation of error amplifiers . 


-107 



Figure A-3d 
Simulation of Ramp. 





Simulation of motor field control circuit 





by following a sequential testing procedure as described below. In 
each step some useful results is obtained. In the first three tests 
described below the motor field flux is held at 1 p.u. to simulate 
generator field control mode. 

Step-1 

We open all but the inner most loop (in our case Gen.51d.Curr. loop) 
and apply a step input as reference. This value should preferably 
be larger than that of a unit step so that one is assured that 
voltage amplitudes are well below saturation level of amplifiers 
during transients in the total simulation. Fig.A-4 depicts the 
field current response for a step size of 1.75 p.u. step input. 

In block diagram (Fig.A-2) this input corresponds to input to field 
mixing amplifier (Block A) and physically corresponds to a step 
error in armature current. This response is satisfactory snd does 
not require any compensation. 

Step-2 

Wo now test the enveloping loop (armature current loop) by applying 
a step input as reference and the procedure of step 1» is repeated. 
The uncompensated response is shown in Fig .A. 5^ We find that the 
armature current response has large overshoot. To improve the 
response a feedback signal through a phase lead compensating 
network (known as phase advancer in industry) is given, 
parameters of compensating block are adjusted to obtain satisfactory 
response for both large and small step reference inputs. 
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Step-5 

Next the perforraance of speed control loop is checked. It was found 
necessary to increase voltage feedback gain to 200 to ensure little 
overshoot and regenerative current. Next this circuit is tested for 
Ramp input (Pig.A-6a,b). 

Step-^A 

To check for field weakening the circuit is connected completely and 
reference signal corresponding to run to top speed is given. The 
compensating parameters in the feedback loop from the tertoinal 
voltage are obtained for no load. It was felt that an additional 
l.oop with phase advance will have to be- included to reduce speed drop 
on impact loading and load throwoff (ref. Pig.A-ya-d). 

■ w 

A test run with I .5 p.u. full load is taken to ensure that during 
short term and accidental overloads the mill does not fail to operate. 

The summary of the overall performance is given in Table A-1 . 
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Tatle A-1 


Geu.Pld.Ckt. 

With forcing function 

Without forcing function 

Time Const. 

0.1875 seconds 

3.44 seconds 


Arni8 ture 

With phase advance 

Without’ phase advance 

Circuit 
Overshoot 
(Step Input) 

Large step 

5.6/« 

17.6^ 

Sma.ll step 

1.59^ 



_ , Term. Vol. Feedback gain 200 60 

Speed loop 

performance Overshoot almost 259^ 

(stop Input ) 

Regenerative current 0.02 p.u. 0,4 p.u. 



'^’l’'ith vreak 
field 

Impact 

Recovery Arm. 
time Current 

(Less than) (p.u.) 

Arm.Curr, Term. S.S.Top 
peak Voltage speed 

(p.u.) io (p.u.) 

Load On 

109^ 

1.5 sec. 

0.98 

1.1 10 

1.9 

Load 

Throwoff 

5CP/0 

5 sec. 

avg.O 

1.8 10 

2.1 



APPENDIX - B 


Ihe Natxge of Stalling Torque dixTing Rolling 

During rolling, the mill drive is required to overcome torque due to 

the metal slippage past the arc of contact a) with rolls and 

deformation of steel tshich is in plastic condition. Ihe mill drive 

8 1 1 

experiences this torque as a stalling - torque ’ which is evident 
from the theoritical torque vs. time charac't'erstics shown in 
Pig.B-1t. The gentle upward slope in the characteristics of 
Pig.B-lb indicates higher torque requirement as rolling proceeds 
primarily because of cooling (radiation loss) of the metal. It is 
assumed throughout this investigation that torque vs. time , 
characteristics are prespecified by a process’ computation from 
slab dimension and furnace data. However, for the purpose of 
present simul-ation the characteristics as shown in Pig.B-1c has 


been used. 




Figure B-1c 
Figure B-1 

The nature of torque during rolling. 




APPENDIX - C 


The Derivation of Bellman's Functional Equation 
Let the state equations for the system he 

X - f(x, u, t) C.1.1 

where x is n x 1 vector defining the state of the system 
u is r X 1 vector defining the forcing functions 

x(0) « x(o), t ^ (O, T) 

It is desired to extremise 

P(x, u, t) dt C,1.2 

./q 

For such problems, Bellman has siiggested a direct computational 
procedure v/hich is a consequence of imbedding principle and optimality 
principle , 

2 

Imbedding Principle * A problem with a fixed initial state and a 
fixed interval of operation may always he viewed as a special ease 
of a more general problem, with a variable initial state and variable 
interval of operation. 

If t' is such that 0 ^t!^ T then the original problem ean be 
considered as a particular case of a more general problem vAiexe 
x(t') defines the initial state and (t', T) defines the interval 
of operation. 

Optimality Principle » Any choice of u ishich is optimal over the 
interval (O, T) is necessarily optimal over any sub-interval (t%t‘') 



where 0 <r'"t ' ^ T . 

\' ^ 

V , 

Thus the optimum value of E is a function of initial state and time, 
o p T 

E(x(o), O) = optimum E = min 1 i P(x, u, t) 

u e U I 0 

i.. 

Then by invoicing optimality principle and imbedding principle, 

we have 

E(x(o) 50) = min ( I . P(x,u,t) dt +E(x(t'), t') : C.I.4 

u U 0^ ^ J 

0 

Let t' = A t, then because E is a continuous function of x and t 
we have 

E(x(tO,t) = E(x( 0),0) + Udi(x(0),0)/^ t)+ 

+ higher order terms .. .. G.1,5* 

Tlio substitution of C.I.5 in 0,1,4 yields the necessary condition 
for optimality 0,1.6 

min I (>^i/;^t) + C^> i/A x) * X + F j = 0 .. 0 . 1.6 

u c" U L ■■ i 

The equation 0 . 1 . 6 is knovn as Bellman's fxmctional equation. 



* Transposition is implied when ever necessary. 






